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1 Introduction 

In understanding the dynamics of the kaon nucleon scattering or photo-kaon production 
in nucleon, it is important to know the hadronic coupling constants involving the kaons. 
Among them, gx\N and Qky,n are the most relevant coupling constants. Phenomenological 
models for determination of these constants from kaon-nucleon scattering and from the kaon 
photo-production, involve many unknown parameters (see for example |l| and references 
therein). Therefore any prediction about these constants is strongly model dependent and 
suffers from large uncertainties. For this reason a quantitative calculation of the gxYN (Y = 
A or E) coupling constants with a tractable and reliable theoretical approach is needed. 

It is widely accepted that QCD is the underlying theory of the strong interactions. In the 
typical hadronic scale the strong coupling constant a s (/x = rrihad) becomes large and QCD 
is nonperturbative. For this reason calculation of Qkyn is related to the nonperturbative 
sector of QCD, and some kind of nonperturbative approach is needed for determination of 
the above-mentioned coupling constants. Among various nonperturbative methods, QCD 
sum rules [Q is a powerful one. This method is based on the short distance OPE of 
vacuum vacuum correlation function in terms of condensates. For the processes involving 
light mesons n, K or p, there is an alternative method to the traditional QCD sum rules, 
namely, light cone QCD sum rules 0. In this approach the expansion of the vacuum-meson 
correlator is performed near the light cone in terms of the meson wave functions. The 
meson wave functions are defined by the matrix elements of non-local composite operators 
sandwiched between the meson and vacuum states and classified by their twists, rather 
than dimensions of the operators, as is the case in the traditional sum rules (more about 
application of light-cone QCD sum rules can be found in [|5|-[12| and references therein). 



In this work we employ light cone QCD sum rules method to extract the coupling 
constants Qkyn- These coupling constants were investigated in framework of the traditional 
QCD sum rules method in [[TJ, |TB[ for the structure ^75, and for the structure 0^75 in 
|14|| . The discrepancy between the results of these works makes it necessary to perform 



independent calculations in determining the coupling constants Qkyn- In the present article 
we restrict ourselves to the consideration of the structure 0^75 whose choice is dictated by 
the following reason. In |T5], |16| it was pointed out that there is coupling scheme dependence 
for the structures 75, j/75, (i.e., dependence on the pseudoscalar or pseudovector forms of the 
effective interaction Lagrangian of pion with hadrons in the phenomenological part, have 
been used), while the structure 0^75 is shown to be independent of any coupling schemes. 
In order to calculate the coupling constants gxYN we start with the following two-point 
function 

IL(p,p 1 ,q)=Jd i xe* c (0\T7iy(x)fj N {0)\K(q)) , (1) 

where p and i]y are the four-momentum of the hyperon and its interpolating current, 
respectively, tjn is the nucleon interpolating current, q is the four-momentum of K meson. 
The interpolating currents for A, £ and A" are |17], 18] 
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where s, u and d are strange, up and down quark fields, respectively and C is the charge 
conjugation operator, a, b and c are the color indices. 



As has already been mentioned, it was pointed out in [14], 0| that a better determination 
of g-nNN can be done by the structure 0^75, since this structure is independent of the 
effective models employed in the phenomeno logical part. This fact motivates us to calculate 
9nyk in this structure. 

Using the Lorentz, parity and charge conjugation invariance, T(p,pi,q) can be repre- 
sented in the following general form 

T(p,Pi,q) = T 1 (p 2 ,p 2 1 ,q 2 )-f 5 + T 2 (p 2 ,p 2 1 ,q 2 ) d-f 5 

+T 3 (p 2 ,pl q 2 ) /> 7 5 + Up 2 ,pl q 2 )a, vl ^q v , (3) 

where q = p — p\, P = (p + pi)/2 . 

On the phenomeno logical part, these different Dirac structures are obtained by saturat- 
ing correlator (1) of the Y and N states 



(Y(p)\K(q)N( Pl )) (0\ VY (x)\Y(p)) (N( Pl )\rj N (0)\0) 



(4) 



(p 2 -m Y ){p\ -m 2 N ) 
The matrix elements in Eq. (4) are defined in the following way 

(0\ VY (x)\Y(p)) = \ Y u{p) , 
(iV(pi)|^(x)|0) = \ N u(Pi) , 
(Y(p)\K(q)N( Pl )) = g KY Nu(p)ls<Pi) ■ (5) 

Substituting Eq. (5) in Eq. (4), we get 

T = — 2 y~ (p 1 + m Y ) 75 (p'l + m N ) + higher resonances , 

(p — m Y ){Px — Tn N ) 



which can be written as 

gKYN^Y^N 
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{m Y m N - ppi) 75 H $75 



^75 - io a pPa(lplb ■ (6) 



2 
Choosing the structure cr Q/ g7 5 as the physical part, we have 

T*» = -i 9KY 2 N ^ Y 2 XN Pa9fi . (7) 

(p 2 -m Y )(p( -my) 

Let us now turn our attention to the theoretical part of the correlator (1). From Eq. (1) 
we immediately get 



T = aJdxe^ x {-4^ fl tS^ 5 (0\u(0) lu aS T C- 1 lfl s(x)\K(q)] 
=F757m^7^757m^7^75 (0 \u(0)^ 5 s{x)\K(q)) 
T757m^7 v 7 5 7p7m^7^75 (0 \u(0)j p j 5 s(x)\K(q)) } , (8) 



where upper (lower) sign corresponds to A (E) case, a = y2/3 (v^2) for A (£). In Eq. 
(8) S is the full light quark propagator containing both perturbative and nonperturbative 
contributions, 
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(9) 



It follows from Eq. (8) that, in order to calculate the correlator function in QCD, the 
matrix elements of the nonlocal operators between the vacuum and kaon states are needed. 
These matrix elements are defined in terms of kaon wave functions, and up to twist four 
these wave functions can be written as || [7J: 



{Q\u(Q)^ 5 s(x)\K(q)) = if w qj due- iu « x [y K {u) + x 2 9l (u)] 
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where 
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and / Dai = / da± da 2 «3 5(1 — o.\ — a 2 — CK3) 



(10) 



Due to the choice of the gauge x^A^(x) = 0, the path ordered gauge factor Ve^ 9s ■> ux ^ ux >) 
has been omitted. The wave function <pk{u) is the leading twist r = 2, gi(u), ip±, ip\\, 0j_, (p\\ 
are twist r = 4, and <p p (u) and <f a {u) are the twist r = 3 wave functions. Using Eqs. (8), 



(9) and (10), for the structure da^ we get 
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where 



G,(u) 



g 2 (v)dv . 



In deriving this equation we omit the terms which are equal to zero after integration over 
x. After Fourier transformation for the theoretical part of the correlator function, we get 



T theor = -i a f K p a qfiUjdu(p K (u) 
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(12) 



According to the general strategy of QCD sum rules, the quantitative prediction for Qkyn 
can be obtained by matching the representations of a correlator (1) in terms of hadronic 
(Eq. (7)) and quark-gluon degrees of freedom (Eq. (12)). Equating Eq. (7) and Eq. (12), 
and performing double Borel transformation for the variables p 2 and p\ in order to suppress 
higher state and continuum contributions, we finally get the following sum rules for gx\N 
and gxnN coupling constants: 
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(13) 
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where the function 



/n(So/M 2 



so/M 2 



E 

fc=0 



(^o/M 2 ) fc 
fc! 



is the factor used to subtract the continuum, which is modeled by the dispersion integral in 
the region s±, S2 > so, so being the continuum threshold (of course the continuum threshold 
for Eq. (13) is different than that for Eq. (14)), 
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and M\ and M\ are the Borel parameters. Since masses of N, A and S are very close to 
each other, we can choose M 2 and Mf to be equal to each other, i.e., M\ = M\ = 2M 2 , 
from which it follows that u$ = 1/2. 

From Eqs. (13) and (14) we see that the coupling constants Qkkn and Qkt,n are de- 
termined by the quark condensate and wave functions (for the structure cr a/ 37 5 ). We can 
deduce from these expressions that, the coupling constant Qkkn is determined mainly by the 
lowest twist (r = 2) wave function ipx(u), and the gxT,N is determined by the twist (r = 4) 
wave function, hence we expect that gx\N > Qkhn- Indeed our numerical calculations 
confirm this expectation, as is presented in the next section. 

2 Numerical analysis 

The principal nonperturbative inputs in the sum rules (13) and (14) are the kaon wave 
functions on the light cone. In [[| a theoretical framework has been developed to study 
these functions. In particular, it has been shown that the wave functions can be expanded 
in terms of the matrix elements of conformal operators which in a leading logarithmic 
approximation do not mix under renormalization. For details, we refer the reader to the 
original literature @, |j|. In our numerical analysis we use the set of wave functions proposed 
in |0. The explicit expressions of wave functions and the values of the various parameters 
are: 
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where the C 2 and G 4 are the Gegenbauer polynomials defined as 
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5(2u-l) 2 ±l 



C , 4 / (2m-1) = — 21(2u-l) 4 -14(2u-l) 2 + l 



(16) 



and a 2 (fi = 0.5 GeV) = 2/3 and a^n = 0.5 GeV) = 0.43. The parameter <5(u) 2 was 
estimated from QCD sum rules to have the value £ 2 (u) = 0.2 GeV 2 fl9| , e(u = 1 GeV) = 
0.5 §. Furthermore we take /* = 0.156 GeV, fjL K (jjt = 1 GeV) = 1 GeV, m 2 = 
0.8 GeV 2 , (gg)U=i Ge y = -(0.243 GeV) 3 , 4 = (m A + 0.5) 2 GeV 2 , s$ = (m E + 0.5) 2 GeV 2 . 
Also remember that all further calculations are performed at u = uq = 1/2. 

Having fixed the input parameter, one must find the range of values of M 2 for which 
the sum rules (13) and (14) are reliable. The lowest possible value of M 2 is determined 
by the requirement that the terms proportional to the highest inverse power of the Borel 
parameters stay reasonable small. The upper bound of M 2 is determined by demanding 
that the continuum contribution is not too large. The interval of M 2 which satisfies both 
conditions is 1 GeV 2 < M 2 < 2 GeV 2 . The dependence of Eqs. (13) and (14) on M 2 is 
depicted in Figs. 1 and 2. From these figures one can directly predict 



gKAN^A^N 



-0.008 ±0.001 , 
-0.0006 ±0.0001 . 



;i7) 



In determining the values of the strong coupling constants gxAN and Qkhn we need 
the residues of hadronic currents, i.e., A at, Aa and As, whose values are obtained from the 
corresponding mass sum rules for the nucleon, A and £ hyperons [{P7| , fLS}] , as follows 
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(19) 

(20) 
(21) 



where 

a = -2tt 2 (qq) , 

b = ^l^o.uGeV 4 , 

7 = M-l = -0.2, 

(ss) 

and the functions f2(x), fo(x) describe subtraction of the continuum contributions, whose 
explicit forms are presented just after Eq. (14). Dividing both sides of Eq. (17) \a\n and 
Eq. (18) by A^A^, whose numerical values are obtained from Eqs. (19), (20) and (21), 
respectively, for gxAN and gxzN coupling constants we get 

\9kan\ = 10 ± 2 

\9kzn\ = 0.75 ±0.15 (22) 

Let us compare our predictions of gKAN and Qkt,n coupling constants with that of 
traditional sum rules results for the structure a^vlhVyflv 01 ■ The resu lts for these quantities 
in framework of the traditional QCD sum rules method are 

\g KAN \ = 2.37 ±0.09, 

\9kxn\ = 0.025 ±0.015. (23) 

When Eqs. (23) and (24) are compared, it is observed that the light cone predictions 
on gKAN and gxY.N are approximately 4 and 30 times larger, respectively, compared to that 
of the traditional QCD sum rules results. 

As an additional remark, it should be noted that the values of the coupling constants 
gKAN and gxT,N obtained in this work differ from that of the exact SU(3) prediction. Using 
de Swarfs convention PC| , SU(3) symmetry predicts 

gKAN = — at(3 - 2av)g-!tNN , (24) 

9KT.N = (2«d - l)g n NN ■ (25) 

Taking a-® = 0.64 [21] from Eqs. (24) and (25) we have 

gKAN 



gKSN 



3.55 , (26) 



while our result for this ratio is \gKAN/9KHN\ — 12. 

Finally we would like to state that, a more detailed analysis for determination of the 
above-mentioned coupling constants from different structures (75, /P75, ^75) is needed. 
Such an analysis may help to understand the source of discrepancy between predictions of 
different structures. 



Figure captions 

Fig. 1 The dependence of Qkan^a^n on the Borel parameter M 2 . 
Fig. 2 The same as in Fig. 1, but for gxzN^z^N- 
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